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The large Nc nuclear potential puzzle arose because three- and higher-meson exchange contri- 
butions to the nucleon-nucleon potential did not automatically yield cancellations that make these 
contributions consistent with the general large Nc scaling rules for the potential. Here it is proposed 
that the resolution to this puzzle is that the scaling rules only apply for energy-independent poten- 
tials while all of the cases with apparent inconsistencies were for energy-dependent potentials. It 
' is shown explicitly how energy-dependent potentials can have radically different large Nc behavior 

Othan an equivalent energy-independent one. One class of three-meson graphs is computed in which 
■ the contribution to the energy-independent potential is consistent with the general large Nc rules 

' even though the energy-dependent potential is not. 

(N : 

■ I. INTRODUCTION 

> : 
(N ; 

J The nature of the nucleon-nucleon interaction is at the heart of nuclear physics. One traditional picture of nuclear 
' interactions at low energies is that they are mediated via meson-exchange. From the QCD perspective one can 
0^ envision the quarks and gluons organizing themselves into hadrons and then the baryons interact amongst themselves 
. via the exchange of mesons. It is not immediately clear how one can test this picture how nucleon-nucleon interactions 
' emerge from QCD since we have no a priori method for deriving these interactions directly from QCD. In ref. |^ 
. it was pointed out that large Nc QCD can provide some insight into the issue. Since the meson-exchange picture 
' of the nucleon-nucleon interaction if valid is justified on rather generic grounds, it should be expected to be hold 
"V* ', at any Nc greater than unity and hence should hold for large Nc- However, as noted in ref. |^ it is by no means 
' Q ■ obvious that the meson-exchange picture is in fact consistent with large Nc counting rules. In particular, there is a 
^ ' threat that while a one meson-exchange description yields nucleon-nucleon interactions which are consistent with the 
P5 i large Nc counting rules § |, multiple -meson exchange graphs yield interactions which are not. A consistent large 
Nc description requires a cancellation of all of these dangerous graphs. It was shown in ref. ^ that for all two- 
• i-H : meson exchange such cancellations do in fact occur, provided that the large Nc scaling rules for the nucleon-nucleon 
interaction are interpreted as applying for a nucleon-nucleon potential for use in a Schodinger or Lippmann-Schwingcr 
H ', equation (as opposed to a kernel in a four dimensional Bethe-Salpeter type equation). This result seems to support 
■ - - 1 the traditional meson-exchange picture for nucleon-nucleon potentials. Clearly, this support would be stronger if the 
cancellations seen for two-meson exchange also happen for general multiple-meson exchanges. 

Unfortunately, in ref. Q it was shown that the extension of the techniques of ref. |^ does not automatically lead to 
the types of cancellations seen for two-meson exchanges. Thus, in the absence of some type of conspiracy leading to 
such cancellations there appears to be no way to justify the meson-exchange picture from large Nc QCD. The result 
is puzzling — why should such cancellations occur for two-meson exchange for multiple channels and at the same time 
fail for three-or-higher meson exchange? The purpose of this note is to resolve this puzzle. As will be argued below, 
there is a reorganization of the analysis that yields precisely the type of conspiracy needed to resolve the puzzle. 

To begin with let us consider the background to the problem. Large Nc QCDj^, |^ has proven to be a powerful 
tool to learn about qualitative and semiquantitative features of hadronic physics. In principle one may hope that 
it will also provide important insights into nuclear physics. The possible implications of large Nc QCD for nuclear 
interactions were already evident in Witten's original paper on baryons in large Nc QCD. Witten pointed out that 
(i) the nucleon-nucleon interaction is characteristically of order Nc, (ii) that nucleon-nucleon scattering observables 
had no smooth limit as Nc — > oo if the nucleon momenta were of order Nc, and (iii) for momenta of order Nc, a 
relativistic time-dependent Hartree formalism is appropriate and has a smooth large Nc limit. In practice, no such 
time-dependent Hartree calculations have been carried out, although recently it was shown what type of observables 
are calculable in principle in this framework and the spin- and isospin-dependence of these observables was deduced. 
The problem of nucleon-nucleon interactions for momenta of order iV^ is of real significance. Kaplan and Manohar[Q 
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following the work of Kaplan and Savage Q have argued that useful information about the nucleon-nucleon potential 
can be extracted in this regime. In particular, they suggest that the nucleon-nucleon potential can be associated 
with quark-line connected diagrams between two color-singlet clusters of Nc quark lines. This is motivated by 
Witten's Hartree analysis. Combining this with the known contracted SU(4) spin- flavor symmetry of two flavor 
QCD||, I^, 11, ^ for the coupling to each cluster, this was used to deduce that the large Nc scaling behavior for 



the various spin and isospin contributions to the potential. The leading scaling behavior is given by 

Vi=j ~ Nc ; Vi^j ^ N-' , (1) 

where the subscript indicates the quantum numbers of the exchange in the t channel. It is worth noting that nucleon- 
nucleon potentials which are fit to scattering data have a pattern which is consistent with eq. (|l|) in the sense that the 
components which are order Nc in eq. (|l|) are characteristically significantly larger than those which are of order 1/Nc 
[Q, |l3j . Recently there have also been attempts to study this regime from the perspective of effective field theory 

The issue of consistency is studied in the following way. First one supposes that there exists some hadronic field 
theory whose masses and couplings scale with Nc according to the standard large Nc rules It was pointed out 
in refs.||, |l^ that the large Nc scaling rules of meson-baryon couplings given by the contracted SU(4) spin-flavor 
symmetry of two flavor QCD|||, ||, |o| yield a one- meson exchange potential and will satisfy eq. (|l|). At two- meson 
exchange the key issue is that both the box-graph and the crossed-box graph are formally of order Nc and hence 
individually cannot be consistent with eq. (|l|). It is straightforward to show, however that the box graph can be 
decomposed into two parts — a contribution arising from the nucleon poles and a contribution arising from the meson 
poles when the graph is evaluated via contour integration. The contribution arising from the nucleon poles can be 
shown to be of exactly the same form as an iterate of the one-meson exchange potential in a Lippmann-Schwinger 
equation. Thus, the nucleon pole contribution to the graph will be picked up when solving the Schrodinger or 
Lippmann-Schwinger equation and must not be included as part of the potential in order to avoid double counting. 
The meson-pole contributions to the box graph are retardation effects. The full contribution to the potential from 
these graphs is the sum of the retardation part of the box-graph with the crossed-box. 

In ref . , it was shown that for all types of two- meson exchange there are cancellations between these two graphs 
so that the sum is consistent with eq. (|l|). These cancellations were highly nontrivial since a number of different 
spin and isospin structures must all cancel. For example, two-pion exchange, for which both the box and crossed-box 
contributions are of order N^, contribute to, among other things, the isoscalar central potential, which is of order 
Nc and requires cancellations of at least relative order 1/Nc and to the isovector central potential, which is of order 
1/Nc and requires cancellations of a least relative order 1/Nc, order to maintain consistency. By explicitly checking 
meson exchanges for all relevant spin and isospin couplings it was seen that all of the "dangerous" contributions that 
were inconsistent with eq. (|l|) canceled to the degree necessary to ensure consistency. This demonstration required 
explicit use of the contracted SU(4) algebra which in turn implied that intermediate A states have to be kept as 
explicit degrees of freedom in the potential model to obtain consistency. 

The analysis of ref. |^ clearly helps justify the meson exchange picture of nucleon-nucleon forces. However, no 
general theorem was proved. Rather all of the relevant cases for two-meson exchange were individually tested. 
Clearly, one's confidence in the generality of the result would increase if a number of examples of three- and higher- 
meson exchange show the same type of cancellations seen in two-meson exchange. An analysis of certain multi-meson 
exchange graphs was done in ref. [|6). The diagrams involved can get quite complicated. Accordingly, it was necessary 
to establish some bookkeeping rules. The basic method used was in many ways analogous to that used in ref. 
First one identifies all two-baryon irreducible Feynman diagrams as contributing to the potential. Next one considers 
the two baryon-reducible graphs and notes that such graphs all contain parts that have two baryon propagators 
between interactions. These propagators are then expressed as the sum of two parts — one where one of the nucleons 
is on-shell, and the remainder. Next, one makes use of the fact the two baryon propagators with one baryon on-shell 
is identical to the propagator in the context of a Lippmann-Schwinger equation (up to relativistic corrections that are 
suppressed by 1/Nc). Thus, these contributions will be included when iterating the Lippmann-Schwinger equation 
the potential is identified as coming from the two-baryon reducible graphs as the full graph minus the contributions 
arising from two propagators between interactions with one baryon on shell. 

This organization of the full problem into a potential and then its iteration via a Lippmann-Schwinger or Schrodinger 
equation has a number of virtues. First, provided the problem is nonrelativistic, the scattering amplitude obtained 
by such a procedure correctly reproduces the sum of all Feynman diagrams for on-shcU scattering. Secondly this 
procedure is well suited to the study of ladder and crossed-ladder graphs including various non-commuting couplings 
since the non-abelian generalization of the eikonal formula of refs. |jl^, |l^ can be implemented for the sum of these 
graphs. The non-abelian generalization of the eikonal formula expresses the sum of all meson and crossed-meson lines 
entering a single baryon in terms of commutators multiplying S function which correspond to an on-shell baryon. In 
the context of summing ladder and crossed-ladder contributions to baryon-baryon scattering it is straightforward to 
see which of these on-shell contributions correspond to iterates of the Lippmann-Schwinger equation. 
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FIG. 1: Three scalar-isoscalar meson exchange diagrams with one seagull type coupling. 



When this organizing principle was implemented for multi-meson exchanges, however, it was found that the cancel- 
lations needed for the potentials to maintain consistency with eq. did not occur. This was seen for two classes of 
graphs that were considered. One class was the sum of ladders and crossed-ladders with non-commutating couplings. 
The most dangerous type was where the meson coupled to baryon in a vector-isovector manner. The issue there 
was the emergence of contributions in which one nucleon was on-shell but which were not iterates of the potential as 
defined above. The generalized eikonial formula||l^, |l^ implies that these yield contributions to the potential given in 
terms of commutators. Although the commutators typically yield a ^/N"^ suppression and may be expected to induce 
a commutator on the other nucleon leg for another factor of the ladder has six rungs or higher, the explicit 

iVc associated with each meson-nucleon vertex overwhelms the suppression and an inconsistency with eq. (|^) is the 
result. A simpler case where an inconsistency can be seen is the case of three scalar-isoscalar meson exchange between 
nucleons where two of the mesons couple to one of the nucleons in a seagull type vertex as shown in fig. (^. Again 
when the part representing an iterate of the potential using the organizing principle discussed above is removed, the 
remaining contribution which contributes to the potential itself does not cancel and is of order N'^, in violation of 
eq. (P. This inconsistency is the large Nc potential puzzle. 

In the remainder of this article a resolution to the puzzle will discussed. In the following section, it will be 
suggested that the heart of the problem lies in the organizing principle discussed above which has the feature that the 
potential obtained have explicit energy dependence. Following this a toy problem will presented to show how energy 
dependence can alter the the counting of a potential. Next a new organizing principle is suggested which yields 
energy independent potentials. In the final section it is shown explicitly that this new procedure, when applied to the 
case of three scalar-isoscalar meson exchange where two of the mesons couple to one of the nucleons in a seagull type 
vertex (of the type in fig. (|l|)), yields a potential consistent with eq. (|l|). 



II. ENERGY DEPENDENCE 

In ref. |^ various possible resolutions to the puzzle were suggested. One possibility was that necessary cancellations 
might happen naturally and generically (without any special conspiracies involving coupling constants of different 
mesons) with a different organization of the problem. However, no plausible reorganization was suggested. The 
purpose of the present paper is to argue that this is the correct resolution of the problem and to provide the needed 
reorganization. It will then be shown explicitly that with this reorganization, the contributions to the potential of the 
three-meson exchange graphs of fig. (0) indeed give rise to the cancellations needed for the consistency with eq. (0). 

The key to this reorganization is the realization that the separation of the contributions into a potential and its 
iterates in a Lippmann-Schwinger equation is not unique if one allows energy-dependent potentials: there are an 
infinite number of ways to distribute the energy dependence between explicit energy dependence of the potential and 
energy dependence arising from iteration which yield identical scattering amplitudes. Of course, one often thinks of 
potentials as being energy independent. However, energy-dependent potentials arise naturally when one suppresses 
explicit inclusion of degrees of freedom and in the present problem the final NN potential suppresses explicit meson 
degrees of freedom that are in the underlying problem. The procedure used in ref. |^ described above to isolate the 
potential from its iterates produces such an energy-dependent potential. Of course, at a fundamental level, there is 
nothing wrong with energy dependent potentials for use in a Schrodinger equation provided they correctly predict 
the physical scattering amplitudes, and the procedure used in ref. should reproduce the scattering amplitude. The 
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notion that the eUmination of en ergy dependence to estabhsh consistency may remind one of the situation in NRQCD 
where field redifintions are used. ||l8| 

Since the correct physical amplitude is obtained, the fact that the procedure used in ref. ([^) produces an energy- 
dependent potential might seem completely innocuous. However, the key point is that the large Nc scaling rules in 
eq. (0) is not for a physical quantity which can be directly measured; rather it is for the potential. Moreover, it is 
possible for an energy-dependent potential to have radically different Nc scaling behavior from an energy-independent 
potential which is physically equivalent. 



III. A TOY PROBLEM 



To see how energy dependence can alter the Nc scaling consider the following simple example: Begin with a simple 
energy-dependent potential operator, 

ViE) ^Vo + ViiE- f/MN) , (2) 

where the tilde is used to distinguish the potential from an equivalent energy-independent one; Vq and Vi have no 
energy dependence (so that all energy dependence is explicitly given), p denotes the relative momentum operator, and 
Mn/2 is the reduced mass. The potential is not explicitly Hermitian but this is not significant for the present purpose, 
which is merely illustrative. This nonhermitian form is used since it is the simplest model which demonstrates the key 
point. The scattering amplitude T is obtained via the Lippmann-Schwinger equation, which as an operator equation 
can be represented as 

T = V + VGlsT = V + VGlsV + VGlsVGlsV + . . . 

with Gls = {E- p^/Mn + if)-^ . (3) 

Now suppose that we have another potential V which is energy independent and has the same on-shell scattering 
amplitude as in eq. By comparing the iterates of eq. (||) for the two potentials, it is straightforward to see that 

V = Vo + ViV(, + ViVi% + ... J2 ^iVo- (4) 

j=0,oo 

The products of the Vi and Vq emerge in V from iterates of the Lippmann-Schwinger equation for V; the {E — fP' /Mm) 
factors accompanying Vi when multiplied by G^s give unity and thus correspond to an uniterated term in the 
Lippmann-Schwinger equation for V . Also, the [E — /Mm) operators annihilates the on-shell external state and 
thereby eliminates many terms in the sum. Thus, we see explicitly that energy dependence reshuffles what goes into 
the potential and what is obtained via iteration. For the present context this is significant in terms of the Nc behavior. 
From the form of eq. (^), it is easy to see that the energy-dependent potential V and the energy- independent potential 
V cannot both generically be of order Nc. If V is generically of order Nc {i.e. both Vq and Vi are order Nc), then we 
see that V contains all powers of Nc- Conversely, if V and Vi are each of order N^ then Vq contains all powers of Nc. 



IV. A NEW ORGANIZING PRINCIPLE 



We see from the preceding example that energy-dependent potentials need not have the same Nc dependence as 
energy-independent ones. The large Nc scaling rules of eq. ([T]) were derived by associating the quark-line connected 
diagrams with "the potential" . This association is somewhat heuristic and it is not immediately obvious whether it is 
supposed to apply for energy-dependent or energy- independent potentials. It is a reasonable hypothesis, however, that 
it applies to energy-independent potentials. Thus, it is plausible that the failure of the analysis of ref. Q to reproduce 
the Nc scaling rules for eq. (0) is because the separation into a potential and its iterates used in the analysis does not 
ensure that the potential is energy independent. To see why the potential so derived can depend on energy, consider 
the algorithm used in the analysis to isolate the potentials. One calculates Feynman diagrams and subtracts off all 
contributions coming from places where one nuclcon out of a pair that separates interacting subdiagrams is on-shell. 
This removes terms that look like Lippmann-Schwinger equation iterates. The issue is simply that the remaining 
subdiagrams may, themselves, be energy dependent. 

To test whether energy- independent potentials from multiple meson exchange are consistent with the rules of eq. (|^) , 
we must first extract the contributions from energy-independent potentials from various meson exchanges graphs. The 
basic algorithm is similar to the one used in ref. |^ with one modification. The potential contributions from two- 
particle irreducible diagrams will be taken as the static limit {i.e. zero energy limit) of the diagram. For two-particle 
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reducible graphs, one subtracts off all contributions coming from places where one nucleon out of a pair that separates 
interacting subdiagrams is on-shell, and where the remaining subdiagrams are replaced by their static limits. 

It is worth noting at this point that at the level of two-meson exchange, the organizing principle of ref. ^ and the 
present one both give results consistent with eq. (|^). The reason for this is quite simple. The two only differ in the 
treatment of the energy dependence of subdiagrams. In the case of two-meson exchange, the subdiagrams are single 
meson exchanges which do have nontrivial energy dependence due to the poles in the meson propagators. However, 
when the external lines are on-shell, the two methods are identical; both methods have the same cancellation between 
retardation effects in the box diagram and the full crossed-box. Of course, one can distinguish between the two-meson 
exchange in the two methods if one or the other external line goes off-shell. However, external line can go off-shell only 
inside larger graphs and the lowest order one can distinguish between the two approaches is at the level of three-meson 
exchange. 



V. TESTING THE Nc SCALING OF ENERGY-INDEPENDENT POTENTIALS 

The modification suggested in the previous section greatly complicates the analysis. For involved cases, such as 
ladders and crossed-ladders with many rungs, it may require a considerable effort to test the consistency of eq. (|l|) 
since the nonabelian generalization of the eikonal formula can no longer be implemented in a straightforward way. 
However, the case of three scalar-isoscalar meson exchange between nucleons where two of the mesons couple to 
one of the nucleons in a seagull type vertex represented in fig. (Q) is tractable. The diagrams in fig. (|l|) include all 
contributions which has the seagull attached to one of the nucleon lines. There is an identical contribution which 
have the seagull attached to the other line. It is useful to group the diagrams into two sets, (a), (b) and (c) as one 
group and (e), (f) and (g) as the other. The cancellations needed to get consistency with eq. (|]) can be shown to 
occur separately in these groups. Now let us consider the contribution from the first three diagrams. The amplitude 
is given by 



iAaM,c = 9lm92m 




X G{p - ki) [G{p + ki)G{p + ki+ k2) + GiP + k2)G{p + ki+ k2) 

+ G{p + k2)G{p + k2 + k:i)]. (5) 

where boldfaced indicates three vectors, and non-boldfaced four vectors, with kj = {ujj,kj); gim and g2m are the 
one-meson and two-meson coupling constants; the initial four-momenta of the two nucleons are p and p; the four- 
momentum transfer is q; Dik) = (fc^ — -I- ie)~^ is the meson propagator; and the fermion propagator is denoted 
by G. For simplicity we will work in the center of mass frame with p = (p^/(2M„),p) and p = (p^/(2M„), —p). In 
fact, the graph as written is ultraviolet divergent. We will assume that it is regulated by some short distance physics 
which acts to cut off the momentum integrals. The details of how this is done is irrelevant for what follows, provided 
the same cutoff procedure is used for all of the graphs. 

Let us note the Nc dependence of the various inputs to this expression: gim ^ Nc^^ , g2m ^ and Mn ~ N^, the 
meson mass in the propagators D is of order N^- Thus, there is an overall prefactor of coming from the coupling 
constants; it is this factor which must somehow be canceled up to relative order N~'^ in order to get consistency with 
eq. (0) . The kinematic regime of interest is intrinsically nonrelativistic since p ~ while M^r ~ ■ In this regime 
it is legitimate to replace the full fermion propagator by 

= r.2/Lr (1 + ^(V^c)) , (6) 

Po ~ P /{2Mn) + 2 

where the 1/Nc corrections come from the nonrelativistic reduction. The propagator contains a recoil correction. 
This will be a 1/Nc correction everywhere except in the vicinity of the propagator's pole. This pole correction is 
relevant only for the piece which looks like a Lippmann-Schwinger equation iterate and, hence, we will drop the recoil 
correction everywhere except for this one contribution. These recoilless propagators depend only on the energy and 
are given by (po + ie)~^ . 
A tedious calculation yields 

3 = 1 \i=l / 3 = 1 \j = l / 
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X D(fci) [S {po - c^i + (p + fci)V(2M^)) G{p + fci)] 

X D{k2)D{h) ^ +0{N,), (7) 

where the 0{Nc) corrections conie from using nonrelativistic propagators and neglecting recoil. The combination 
in square brackets can easily be seen to be the delta function fixing cui times the Lippmann-Schwinger propagator, 

Gls{P,E) = E-pyMj,+te - 

Note that the delta function restricts tui to being of order for these kinematics. The expression in square 
brackets is sensitive to the infrared kinematics and depends on the fact that oJi ^ 0. However, for all parts of the 
expression, setting lo\ to zero will only induce an error of relative size (1/iVc). Substituting zero for lo\ in all parts of 
eq. except for the factor in square brackets gives 



^ GL5(fcl, i?) V2[q - fci) + 0(iV,) 



k 



3 3 / 3 \ / 3 



1 (2^) 



J=2 \j=2 / \j=2 



X D{k2)D{k3)^— . (8) 

UJ2 + le 

Note that Vi is a static one-meson exchange potential while V2 is the static potential for the exchange of two mesons 
with one seagull vertex. Thus, the form of eq. is precisely part of an iterate of the Lippmann-Schwinger equation, 
with Vi as the first iteration and V2 the second. The other ordering comes from the graphs (e), (f) and (g). Since 
the part of iAa.b,c which is of order is a static potential iterate, one concludes the contributions of these graphs 
to the potential is necessarily of order A^^ or less and, hence, is consistent with the Nc scaling rules of eq. (|^). 

Note that had one followed the organization of ref. |^ there would have been a remaining contribution to the 
potential of order . From this we see an explicit example where the energy-dependent potentials extracted from a 
set of Feynman graphs associated with multiple-meson exchange using the algorithm of ref. |^ are inconsistent with 
the large Nc counting rules of refs. 0, |^ while energy-independent ones are consistent. Thus, we see the large Nc 
nuclear potential puzzle has been resolved for this set of diagrams. Although we have not explicitly calculated other 
multiple-meson exchange graphs such as the sum of ladders and crossed-ladders considered in ref. Q due to their 
complexity, it is highly plausible that consistency would be obtained for the energy-dependent potentials for such 
cases. 

In summary, it has been shown that energy-dependent potentials can have different Nc scaling behavior than 
energy-independent potentials. It has been argued that the Nc scaling rules of eq. (|^) should apply only to energy- 
independent potentials. With these energy-dependent potentials one expects consistency. From this perspective, the 
inconsistency of the three- and higher-meson exchange potentials in ref. ^ with eq. (|l|) can be understood as arising 
from the fact that the analysis used in that work yielded energy-dependent potentials. This interpretation is highly 
plausible given the explicit demonstration of consistency for the case of the energy-independent potential associated 
with three-scalar-meson exchange of the type seen in fig. (|^) despite the fact that the energy-dependent potential was 
obtained using the methods of ref. |^ 

This work is supported by the U.S. Department of Energy grant DE-FG02-93ER-40762. The author gratefully 
acknowledges useful discussions with A. Beltisky, B. Gelman, D. Kaplan, A. Manohar and D. Riska. 
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